3.2.53 \(\int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx\) [153]

3.2.53.1 Optimal result
3.2.53.2 Mathematica [C] (verified)
3.2.53.3 Rubi [A] (verified)
3.2.53.4 Maple [A] (verified)
3.2.53.5 Fricas [F]
3.2.53.6 Sympy [F]
3.2.53.7 Maxima [F]
3.2.53.8 Giac [F]
3.2.53.9 Mupad [F(-1)]

3.2.53.1 Optimal result

Integrand size = 19, antiderivative size = 261 \[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=-\frac {b c}{2 d x}+\frac {b c^2 \text {arctanh}(c x)}{2 d}-\frac {a+b \text {arctanh}(c x)}{2 d x^2}+\frac {e (a+b \text {arctanh}(c x))}{d^2 x}-\frac {b c e \log (x)}{d^2}+\frac {a e^2 \log (x)}{d^3}+\frac {e^2 (a+b \text {arctanh}(c x)) \log \left (\frac {2}{1+c x}\right )}{d^3}-\frac {e^2 (a+b \text {arctanh}(c x)) \log \left (\frac {2 c (d+e x)}{(c d+e) (1+c x)}\right )}{d^3}+\frac {b c e \log \left (1-c^2 x^2\right )}{2 d^2}-\frac {b e^2 \operatorname {PolyLog}(2,-c x)}{2 d^3}+\frac {b e^2 \operatorname {PolyLog}(2,c x)}{2 d^3}-\frac {b e^2 \operatorname {PolyLog}\left (2,1-\frac {2}{1+c x}\right )}{2 d^3}+\frac {b e^2 \operatorname {PolyLog}\left (2,1-\frac {2 c (d+e x)}{(c d+e) (1+c x)}\right )}{2 d^3} \]

output
-1/2*b*c/d/x+1/2*b*c^2*arctanh(c*x)/d+1/2*(-a-b*arctanh(c*x))/d/x^2+e*(a+b 
*arctanh(c*x))/d^2/x-b*c*e*ln(x)/d^2+a*e^2*ln(x)/d^3+e^2*(a+b*arctanh(c*x) 
)*ln(2/(c*x+1))/d^3-e^2*(a+b*arctanh(c*x))*ln(2*c*(e*x+d)/(c*d+e)/(c*x+1)) 
/d^3+1/2*b*c*e*ln(-c^2*x^2+1)/d^2-1/2*b*e^2*polylog(2,-c*x)/d^3+1/2*b*e^2* 
polylog(2,c*x)/d^3-1/2*b*e^2*polylog(2,1-2/(c*x+1))/d^3+1/2*b*e^2*polylog( 
2,1-2*c*(e*x+d)/(c*d+e)/(c*x+1))/d^3
 
3.2.53.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 5.93 (sec) , antiderivative size = 375, normalized size of antiderivative = 1.44 \[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=-\frac {\frac {a d^3}{x^2}-\frac {2 a d^2 e}{x}-\frac {b d^3 \left (-1+c^2 x^2\right ) \text {arctanh}(c x)}{x^2}+\frac {b e^3 \text {arctanh}(c x)^2}{c}+\frac {b d^2 (c d-2 e \text {arctanh}(c x))}{x}-2 a d e^2 \log (x)+2 a d e^2 \log (d+e x)+b c d^2 e \left (2 \log (c x)-\log \left (1-c^2 x^2\right )\right )-b d e^2 \left (\text {arctanh}(c x) \left (\text {arctanh}(c x)+2 \log \left (1-e^{-2 \text {arctanh}(c x)}\right )\right )-\operatorname {PolyLog}\left (2,e^{-2 \text {arctanh}(c x)}\right )\right )+\frac {b e^2 \left (-\sqrt {1-\frac {c^2 d^2}{e^2}} e e^{-\text {arctanh}\left (\frac {c d}{e}\right )} \text {arctanh}(c x)^2+c d \text {arctanh}(c x) \left (i \pi +2 \text {arctanh}\left (\frac {c d}{e}\right )+2 \log \left (1-e^{-2 \left (\text {arctanh}\left (\frac {c d}{e}\right )+\text {arctanh}(c x)\right )}\right )\right )-\frac {1}{2} i c d \pi \left (2 \log \left (1+e^{2 \text {arctanh}(c x)}\right )+\log \left (1-c^2 x^2\right )\right )+2 c d \text {arctanh}\left (\frac {c d}{e}\right ) \left (\log \left (1-e^{-2 \left (\text {arctanh}\left (\frac {c d}{e}\right )+\text {arctanh}(c x)\right )}\right )-\log \left (i \sinh \left (\text {arctanh}\left (\frac {c d}{e}\right )+\text {arctanh}(c x)\right )\right )\right )-c d \operatorname {PolyLog}\left (2,e^{-2 \left (\text {arctanh}\left (\frac {c d}{e}\right )+\text {arctanh}(c x)\right )}\right )\right )}{c}}{2 d^4} \]

input
Integrate[(a + b*ArcTanh[c*x])/(x^3*(d + e*x)),x]
 
output
-1/2*((a*d^3)/x^2 - (2*a*d^2*e)/x - (b*d^3*(-1 + c^2*x^2)*ArcTanh[c*x])/x^ 
2 + (b*e^3*ArcTanh[c*x]^2)/c + (b*d^2*(c*d - 2*e*ArcTanh[c*x]))/x - 2*a*d* 
e^2*Log[x] + 2*a*d*e^2*Log[d + e*x] + b*c*d^2*e*(2*Log[c*x] - Log[1 - c^2* 
x^2]) - b*d*e^2*(ArcTanh[c*x]*(ArcTanh[c*x] + 2*Log[1 - E^(-2*ArcTanh[c*x] 
)]) - PolyLog[2, E^(-2*ArcTanh[c*x])]) + (b*e^2*(-((Sqrt[1 - (c^2*d^2)/e^2 
]*e*ArcTanh[c*x]^2)/E^ArcTanh[(c*d)/e]) + c*d*ArcTanh[c*x]*(I*Pi + 2*ArcTa 
nh[(c*d)/e] + 2*Log[1 - E^(-2*(ArcTanh[(c*d)/e] + ArcTanh[c*x]))]) - (I/2) 
*c*d*Pi*(2*Log[1 + E^(2*ArcTanh[c*x])] + Log[1 - c^2*x^2]) + 2*c*d*ArcTanh 
[(c*d)/e]*(Log[1 - E^(-2*(ArcTanh[(c*d)/e] + ArcTanh[c*x]))] - Log[I*Sinh[ 
ArcTanh[(c*d)/e] + ArcTanh[c*x]]]) - c*d*PolyLog[2, E^(-2*(ArcTanh[(c*d)/e 
] + ArcTanh[c*x]))]))/c)/d^4
 
3.2.53.3 Rubi [A] (verified)

Time = 0.56 (sec) , antiderivative size = 261, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.105, Rules used = {6502, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx\)

\(\Big \downarrow \) 6502

\(\displaystyle \int \left (-\frac {e^3 (a+b \text {arctanh}(c x))}{d^3 (d+e x)}+\frac {e^2 (a+b \text {arctanh}(c x))}{d^3 x}-\frac {e (a+b \text {arctanh}(c x))}{d^2 x^2}+\frac {a+b \text {arctanh}(c x)}{d x^3}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {e^2 \log \left (\frac {2}{c x+1}\right ) (a+b \text {arctanh}(c x))}{d^3}-\frac {e^2 (a+b \text {arctanh}(c x)) \log \left (\frac {2 c (d+e x)}{(c x+1) (c d+e)}\right )}{d^3}+\frac {e (a+b \text {arctanh}(c x))}{d^2 x}-\frac {a+b \text {arctanh}(c x)}{2 d x^2}+\frac {a e^2 \log (x)}{d^3}+\frac {b c^2 \text {arctanh}(c x)}{2 d}+\frac {b c e \log \left (1-c^2 x^2\right )}{2 d^2}-\frac {b e^2 \operatorname {PolyLog}(2,-c x)}{2 d^3}+\frac {b e^2 \operatorname {PolyLog}(2,c x)}{2 d^3}-\frac {b e^2 \operatorname {PolyLog}\left (2,1-\frac {2}{c x+1}\right )}{2 d^3}+\frac {b e^2 \operatorname {PolyLog}\left (2,1-\frac {2 c (d+e x)}{(c d+e) (c x+1)}\right )}{2 d^3}-\frac {b c e \log (x)}{d^2}-\frac {b c}{2 d x}\)

input
Int[(a + b*ArcTanh[c*x])/(x^3*(d + e*x)),x]
 
output
-1/2*(b*c)/(d*x) + (b*c^2*ArcTanh[c*x])/(2*d) - (a + b*ArcTanh[c*x])/(2*d* 
x^2) + (e*(a + b*ArcTanh[c*x]))/(d^2*x) - (b*c*e*Log[x])/d^2 + (a*e^2*Log[ 
x])/d^3 + (e^2*(a + b*ArcTanh[c*x])*Log[2/(1 + c*x)])/d^3 - (e^2*(a + b*Ar 
cTanh[c*x])*Log[(2*c*(d + e*x))/((c*d + e)*(1 + c*x))])/d^3 + (b*c*e*Log[1 
 - c^2*x^2])/(2*d^2) - (b*e^2*PolyLog[2, -(c*x)])/(2*d^3) + (b*e^2*PolyLog 
[2, c*x])/(2*d^3) - (b*e^2*PolyLog[2, 1 - 2/(1 + c*x)])/(2*d^3) + (b*e^2*P 
olyLog[2, 1 - (2*c*(d + e*x))/((c*d + e)*(1 + c*x))])/(2*d^3)
 

3.2.53.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 6502
Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_.)*((f_.)*(x_))^(m_.)*((d_) + (e 
_.)*(x_))^(q_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*ArcTanh[c*x])^p, ( 
f*x)^m*(d + e*x)^q, x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && IGtQ[p, 0] 
 && IntegerQ[q] && (GtQ[q, 0] || NeQ[a, 0] || IntegerQ[m])
 
3.2.53.4 Maple [A] (verified)

Time = 1.14 (sec) , antiderivative size = 322, normalized size of antiderivative = 1.23

method result size
parts \(a \left (-\frac {e^{2} \ln \left (e x +d \right )}{d^{3}}-\frac {1}{2 d \,x^{2}}+\frac {e^{2} \ln \left (x \right )}{d^{3}}+\frac {e}{d^{2} x}\right )+b \,c^{2} \left (-\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (e c x +c d \right )}{c^{2} d^{3}}-\frac {\operatorname {arctanh}\left (c x \right )}{2 d \,c^{2} x^{2}}+\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (c x \right )}{c^{2} d^{3}}+\frac {\operatorname {arctanh}\left (c x \right ) e}{c^{2} d^{2} x}-\frac {c \left (\frac {2 e \left (-\frac {e \left (\operatorname {dilog}\left (\frac {e c x +e}{-c d +e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x +e}{-c d +e}\right )\right )}{2}+\frac {e \left (\operatorname {dilog}\left (\frac {e c x -e}{-c d -e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x -e}{-c d -e}\right )\right )}{2}\right )}{d^{3} c^{3}}-\frac {2 e^{2} \left (-\frac {\operatorname {dilog}\left (c x +1\right )}{2}-\frac {\ln \left (c x \right ) \ln \left (c x +1\right )}{2}-\frac {\operatorname {dilog}\left (c x \right )}{2}\right )}{d^{3} c^{3}}-\frac {\left (\frac {c d}{2}+e \right ) \ln \left (c x +1\right )+\left (-\frac {c d}{2}+e \right ) \ln \left (c x -1\right )-\frac {d}{x}-2 e \ln \left (c x \right )}{c^{2} d^{2}}\right )}{2}\right )\) \(322\)
derivativedivides \(c^{2} \left (-\frac {a}{2 d \,c^{2} x^{2}}+\frac {a \,e^{2} \ln \left (c x \right )}{c^{2} d^{3}}+\frac {a e}{c^{2} d^{2} x}-\frac {a \,e^{2} \ln \left (e c x +c d \right )}{c^{2} d^{3}}+b c \left (-\frac {\operatorname {arctanh}\left (c x \right )}{2 d \,c^{3} x^{2}}+\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (c x \right )}{d^{3} c^{3}}+\frac {\operatorname {arctanh}\left (c x \right ) e}{d^{2} c^{3} x}-\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (e c x +c d \right )}{d^{3} c^{3}}-\frac {\left (-\frac {c d}{2}-e \right ) \ln \left (c x +1\right )+\left (\frac {c d}{2}-e \right ) \ln \left (c x -1\right )+2 e \ln \left (c x \right )+\frac {d}{x}}{2 c^{2} d^{2}}+\frac {e^{2} \left (-\frac {\operatorname {dilog}\left (c x +1\right )}{2}-\frac {\ln \left (c x \right ) \ln \left (c x +1\right )}{2}-\frac {\operatorname {dilog}\left (c x \right )}{2}\right )}{d^{3} c^{3}}-\frac {e \left (-\frac {e \left (\operatorname {dilog}\left (\frac {e c x +e}{-c d +e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x +e}{-c d +e}\right )\right )}{2}+\frac {e \left (\operatorname {dilog}\left (\frac {e c x -e}{-c d -e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x -e}{-c d -e}\right )\right )}{2}\right )}{d^{3} c^{3}}\right )\right )\) \(340\)
default \(c^{2} \left (-\frac {a}{2 d \,c^{2} x^{2}}+\frac {a \,e^{2} \ln \left (c x \right )}{c^{2} d^{3}}+\frac {a e}{c^{2} d^{2} x}-\frac {a \,e^{2} \ln \left (e c x +c d \right )}{c^{2} d^{3}}+b c \left (-\frac {\operatorname {arctanh}\left (c x \right )}{2 d \,c^{3} x^{2}}+\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (c x \right )}{d^{3} c^{3}}+\frac {\operatorname {arctanh}\left (c x \right ) e}{d^{2} c^{3} x}-\frac {\operatorname {arctanh}\left (c x \right ) e^{2} \ln \left (e c x +c d \right )}{d^{3} c^{3}}-\frac {\left (-\frac {c d}{2}-e \right ) \ln \left (c x +1\right )+\left (\frac {c d}{2}-e \right ) \ln \left (c x -1\right )+2 e \ln \left (c x \right )+\frac {d}{x}}{2 c^{2} d^{2}}+\frac {e^{2} \left (-\frac {\operatorname {dilog}\left (c x +1\right )}{2}-\frac {\ln \left (c x \right ) \ln \left (c x +1\right )}{2}-\frac {\operatorname {dilog}\left (c x \right )}{2}\right )}{d^{3} c^{3}}-\frac {e \left (-\frac {e \left (\operatorname {dilog}\left (\frac {e c x +e}{-c d +e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x +e}{-c d +e}\right )\right )}{2}+\frac {e \left (\operatorname {dilog}\left (\frac {e c x -e}{-c d -e}\right )+\ln \left (e c x +c d \right ) \ln \left (\frac {e c x -e}{-c d -e}\right )\right )}{2}\right )}{d^{3} c^{3}}\right )\right )\) \(340\)
risch \(-\frac {b c}{2 d x}+\frac {c^{2} b \ln \left (-c x \right )}{4 d}-\frac {c^{2} b \ln \left (-c x +1\right )}{4 d}+\frac {b \ln \left (-c x +1\right )}{4 d \,x^{2}}-\frac {c b e \ln \left (-c x \right )}{2 d^{2}}+\frac {c b e \ln \left (-c x +1\right )}{2 d^{2}}-\frac {b e \ln \left (-c x +1\right )}{2 d^{2} x}+\frac {b \,e^{2} \operatorname {dilog}\left (-c x +1\right )}{2 d^{3}}+\frac {b \,e^{2} \operatorname {dilog}\left (\frac {e \left (-c x +1\right )-c d -e}{-c d -e}\right )}{2 d^{3}}+\frac {b \,e^{2} \ln \left (-c x +1\right ) \ln \left (\frac {e \left (-c x +1\right )-c d -e}{-c d -e}\right )}{2 d^{3}}-\frac {a}{2 d \,x^{2}}+\frac {a e}{d^{2} x}+\frac {a \,e^{2} \ln \left (-c x \right )}{d^{3}}-\frac {a \,e^{2} \ln \left (e \left (-c x +1\right )-c d -e \right )}{d^{3}}-\frac {b \,c^{2} \ln \left (c x \right )}{4 d}+\frac {b \,c^{2} \ln \left (c x +1\right )}{4 d}-\frac {b \ln \left (c x +1\right )}{4 d \,x^{2}}-\frac {b c e \ln \left (c x \right )}{2 d^{2}}+\frac {b c e \ln \left (c x +1\right )}{2 d^{2}}+\frac {b e \ln \left (c x +1\right )}{2 d^{2} x}-\frac {b \,e^{2} \operatorname {dilog}\left (c x +1\right )}{2 d^{3}}-\frac {b \,e^{2} \operatorname {dilog}\left (\frac {e \left (c x +1\right )+c d -e}{c d -e}\right )}{2 d^{3}}-\frac {b \,e^{2} \ln \left (c x +1\right ) \ln \left (\frac {e \left (c x +1\right )+c d -e}{c d -e}\right )}{2 d^{3}}\) \(431\)

input
int((a+b*arctanh(c*x))/x^3/(e*x+d),x,method=_RETURNVERBOSE)
 
output
a*(-e^2/d^3*ln(e*x+d)-1/2/d/x^2+e^2/d^3*ln(x)+e/d^2/x)+b*c^2*(-1/c^2*arcta 
nh(c*x)*e^2/d^3*ln(c*e*x+c*d)-1/2*arctanh(c*x)/d/c^2/x^2+1/c^2*arctanh(c*x 
)*e^2/d^3*ln(c*x)+1/c^2*arctanh(c*x)*e/d^2/x-1/2*c*(2/d^3/c^3*e*(-1/2*e*(d 
ilog((c*e*x+e)/(-c*d+e))+ln(c*e*x+c*d)*ln((c*e*x+e)/(-c*d+e)))+1/2*e*(dilo 
g((c*e*x-e)/(-c*d-e))+ln(c*e*x+c*d)*ln((c*e*x-e)/(-c*d-e))))-2/d^3/c^3*e^2 
*(-1/2*dilog(c*x+1)-1/2*ln(c*x)*ln(c*x+1)-1/2*dilog(c*x))-1/c^2/d^2*((1/2* 
c*d+e)*ln(c*x+1)+(-1/2*c*d+e)*ln(c*x-1)-d/x-2*e*ln(c*x))))
 
3.2.53.5 Fricas [F]

\[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=\int { \frac {b \operatorname {artanh}\left (c x\right ) + a}{{\left (e x + d\right )} x^{3}} \,d x } \]

input
integrate((a+b*arctanh(c*x))/x^3/(e*x+d),x, algorithm="fricas")
 
output
integral((b*arctanh(c*x) + a)/(e*x^4 + d*x^3), x)
 
3.2.53.6 Sympy [F]

\[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=\int \frac {a + b \operatorname {atanh}{\left (c x \right )}}{x^{3} \left (d + e x\right )}\, dx \]

input
integrate((a+b*atanh(c*x))/x**3/(e*x+d),x)
 
output
Integral((a + b*atanh(c*x))/(x**3*(d + e*x)), x)
 
3.2.53.7 Maxima [F]

\[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=\int { \frac {b \operatorname {artanh}\left (c x\right ) + a}{{\left (e x + d\right )} x^{3}} \,d x } \]

input
integrate((a+b*arctanh(c*x))/x^3/(e*x+d),x, algorithm="maxima")
 
output
-1/2*a*(2*e^2*log(e*x + d)/d^3 - 2*e^2*log(x)/d^3 - (2*e*x - d)/(d^2*x^2)) 
 + 1/2*b*integrate((log(c*x + 1) - log(-c*x + 1))/(e*x^4 + d*x^3), x)
 
3.2.53.8 Giac [F]

\[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=\int { \frac {b \operatorname {artanh}\left (c x\right ) + a}{{\left (e x + d\right )} x^{3}} \,d x } \]

input
integrate((a+b*arctanh(c*x))/x^3/(e*x+d),x, algorithm="giac")
 
output
integrate((b*arctanh(c*x) + a)/((e*x + d)*x^3), x)
 
3.2.53.9 Mupad [F(-1)]

Timed out. \[ \int \frac {a+b \text {arctanh}(c x)}{x^3 (d+e x)} \, dx=\int \frac {a+b\,\mathrm {atanh}\left (c\,x\right )}{x^3\,\left (d+e\,x\right )} \,d x \]

input
int((a + b*atanh(c*x))/(x^3*(d + e*x)),x)
 
output
int((a + b*atanh(c*x))/(x^3*(d + e*x)), x)